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I. INTRODUCTION 

Time-dependent signals reveal much more about their sources than do static ones. Thus, the period of a pendulum 
' gives us information about the ratio of the inertia and the restoring force. Frequently, a Fourier spectral analysis of 
a time-dependent signal displays steady peaks that are indicative of several independent physical frequencies. Such 
^ ' information is of course extremely useful, not only in helping us model the phenomenon, but also in narrowing down 
>—i ' the model parameters. The question naturally arises whether we might also be able to extract similar constraints 
Qv^ from unsteady temporal signals. 

] All real life systems are very high dimensional (they have a very large number of 'degrees of freedom'), such as well 
T— I . developed turbulence Often though, the observed behavior is dominated by a few degrees of freedom. In that case 
^ ' we may represent these most important of these with a low dimensional dynamics and hide our ignorance about the 
, other ones by treating them as a stochastic component. In some cases we even get lucky and be able to ignore the 
stochastic component all together. This lecture concerns itself with this type of situation. 

Of course, in real systems one is further plagued by extrinsic physical and observational noise. How large a level of 
noise is admissible depends on many factors, not the least of which is the nature of the dynamical system itself. 
— . ^ If the dynamics is stationary (autonomous system) and the signal is not multi-periodic, then it is must be chaotic. 
. Sometimes the time-dependence is much slower than the variability of the source, and for practical purposes we may 
consider the system to be stationary. We refer the student to the many excellent introductions to chaos ||^, ||], j^] 
that abound in the literature. 

'"^ ' Another difficulty of principle arises. In variable star Astronomy, for example, we generally have access only to the 
' [ magnitude of the star (The magnitude is defined as the log of the luminosity m — — 2.5LogL). Can we ever hope to 
extract information about the D-dim dynamics from such a scalar time-series. The answer is perhaps astonishingly, 
pH ■ but fortunately, affirmative. 

. Although the topic is relatively new there are good introductory books and articles available on the problem of 
^ ' extracting a dynamics from a time-series, for example [Q 
• I— I . 

H II. THE UNDERLYING DYNAMICS 

Suppose a physical system can be described by an autonomous system of ordinary differential equations (ODEs) 
with d physical variables that make up the vector Y. 

dV/dt ^G{Y), (2.1) 

(A more complicated higher order system of ODEs can always be cast into this canonical form). We call this system 
a d-dimensional dynamics. With given initial conditions, the solution Y of these ODEs, i.e. the d-dimensional vector 
of the physical phase space that is referred to as the trajectory in physical phase-space. 

In most real situations we are unable to measure Y(i). In fact, we usually measure only one quantity, say m{t). 
(In the case of variable stars this would naturally be the luminosity or perhaps the radial velocity.) This observable 
is generally a complicated and unknown function of m(Y) of the physical state vector. 

We assume that the discrete measurement times {t„} are equally spaced. When this is not the case we have to 
interpolate which unfortunately introduces additional noise. 

The connection between Y(t) and {m„ = m(t„)} will be made through what is termed a phase-space reconstruction. 
From {rrin} we construct the ds-uplets of vectors 

X(i„) = X" = {m(<„), m(i„ ~ A), m(i„ - 2A), . . . m(t„ - {de - 1)A)} (2.2) 
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where the delay A is an integer times the time-intervaL Both the physical space and the reconstruction space are 
Euclidian (d and (Ie are integer). 

We assume now that these vectors themselves are generated by a d^j-dim dynamics 

dx/dt = g(x), (2.3) 

or equivalently by an d^-dim map 

^ F[x"], (2.4) 



It does not matter whether we use the description 2.3 or 2.4 - the two are clearly related with {X" = x(t„)}. 
It is natural to assume that the observable m is a smooth function of the physical phase space variable Y, viz 
{m{tn)} — {m(Y{tn))}- Thus, for example, we assume that the luminosity of a star is a smooth function of the 
dynamical variables, such as radii, velocities and temperatures, for example. 

A theorem now states that there is a one-to-one (differentiable) relation between the vector function Y{t) in the 
d-dim physical phase-space and the vector function x(t) n the dE-dhn reconstruction space, provided that is 
sufficiently large so that the reconstructed trajectory is devoid of intersections and cusps. When the latter condition 
is satisfied one says that the reconstruction space is an embedding space, and no additional benefit can be gained 
from further increasing the dimension ds (A circle can generically be embedded in 3-D; putting it in 4-D or 5-D does 
not change the circle). The theorem furthermore states that ds need not be larger than 2 times the fractal dimension 
of the attractor, rounded off to the next integer value 0. In practical time-series, that one need not go to this bound 
and often ds — d. In the following we will loosely refer to the dimension of the reconstruction space as d^;. Another 
obvious, but important constraint is that the function m(Y) preserves all the necessary information about all the 
degrees of freedom of the dynamics. 

This theorem is extremely powerful because some properties are preserved in the embedding, such as the Lyapunov 
exponents and the fractal dimension [Q 1^ . Once the map has been obtained we can thus infer quantifiable properties 
of an unknown underlying dynamics from the observations of a single variable. This is important because often 
an observational data set is too short or too noisy to allow a direct estimation of these quantities. Of the fractal 
dimensions the Lyapunov dimension di is easiest to calculate, and it is defined in terms of the Lyapunov exponents 
{Xi} as 

where the {Xi} are in ordered in decreasing size and K is the largest integer for which the sum is positive. 

The minimum embedding dimension of course sets an upper bound on the physical dimension d. A lower limit is 
obtained from the obvious fact that a fractal attractor has to be embedded in an integer dimensional space. Thus 

dL<K+l<d^<dE (2.6) 
When we are lucky the two limits coincide and wc can obtain d. 



It is our aim to find a global representation , as opposed to local representations, of the flow, i.e. a single function 
F that evolves the trajectory according to eq. 2.4 in all of phase-space. 

We have no a priori knowledge about how to choose the global map F (or the flow g) . It is however natural to try 
a polynomial expansion. This is particularly appropriate when the dynamics is only weakly nonlinear, as seems to be 
the case for example for stellar pulsations Thus we set 

F(X) = ^CkP,(X), (2.7) 

fe=i 

where the polynomials Pj (X) are constructed to be orthogonal on the data set |^ , , and p is the expansion 
order (typically 4 or 5). 

An equivalent approach is to treat the polynomial expansion as a least squares problem, i.e. one minimizes the 
sum S 
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S = Y,Ec. (2.8) 

Q = l 

E,. = J2{x:+'-F^{X-))' (2.9) 

F(X) = ^B'=Alfc(X) (2.10) 

fc 

with respect to the vector coefficients B'^, where A^k(X) represents all the monomials up to order p that can be 
formed with the components of the vectors X. The sum over n runs over the whole data set. In the following 
we will also use one of the components Ea do measure the error. 

Instead of a standard QR least squares algorithm it is better to use a singular value decomposition (SVD) approach 
to effect the minimization and to obtain the coefficients of the polynomial map ||l^. The SVD algorithm has several 
advantages [ p^ . First, it is a very stable numerical method and it is more than ten times faster than the Gram- 
Schmidt procedure. Second, it allows one to introduce a large number of unknown parameters in the fit (even more 
than there are data points) . This is important because the number of parameters goes up very fast with polynomial 
order p, as (C^"'"'''') and other methods such as QR become singular. With SVD the redundant combinations give 
rise to eigenvalues that fall below the machine error or a self-imposed cut-off tOc value and they are automatically 
eliminated. 

III. A TEST CASE: THE ROSSLER OSCILLATOR 

The system of 3 equations known as the Rossler oscillator [|j is given by 

dxi 

-jT = -^2- X3 (3.1) 
dt 



xi + ax2 (3.2) 



dX2 

~dt 

= b+ {xi- c)x3 (3.3) 

The dimension of the 'physical' phase-space of the Rossler system is thus d—2>. The Rossler oscillator has been widely 
studied § 0. In some regions of parameter space (a, 5, c) the trajectory fills what is commonly a Rossler band (It 
is not really a band but is made out of infinitely fine layers, somewhat like 'pate feuilletee'). 

The trajectory of this system is chaotic {i.e. one of the Lyapunov exponents is positive and nearby trajectories are 
divergent) and the attractor is strange (it has a fractal dimension) Since strange and chaotic behavior seem to 
occur simultaneously for the systems of interest to us we shall simply refer to the attractor as chaotic. 

A. The test-data sets 

For our tests we have chosen the values a=0.2, 6=0.2 and c=4.8 in eq. 3^ for which the solution is chaotic and 
the attractor is in the shape of a band, the so called 'Rossler band' |^|. In the following, when we talk about 
the Rossler band or attractor, we refer to the Rossler system with these specific parameter values. Other parameter 
values give attractors with different properties. For further details of all these tests we refer to Serre et al. [ p2| . 

On top in Fig. 1 we show the temporal segment of xi{t) that we use as test input data in the nonlinear analysis. 
(The data sequence has been obtained by integrating the Rossler system with a fourth order Runge-Kutta algorithm 
with a fine time-step of 0.001.) Since we want to test our method we have specifically chosen a relatively short 
test-data set of 4000 points with 5t chosen so as to give of the order of 60 equally spaced points per 'cycle'. 

We stress that this segment of xi data which is sampled at regular time-intervals, is all the information that we 
allow ourselves in the global fiow reconstruction, i.e. we do not use any information about the other two variables^ 
X2{t) and xy,{t). 

It has been shown that some nonlinear methods can erroneously indicate low dimensional behavior for stochastic 
data ( |l^ ijl^ also cf. the lectures by Provenzale and by Smith in this Volume). It is therefore important that our 
analysis be able to clearly distinguish between the first, deterministic, and the second, stochastic data sets. 

We have thus also created a second, stochastic test-data set as follows iQ. We first 'pre- whiten' the Rossler data set 
by subtracting the dominant frequency and its 8 harmonics, all with their respective Fourier amplitudes and phases 
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- we thus split the signal into to a periodic and an irregular part. The Fourier spectrum of the irregular component 
has a broad band structure, and we compute the envelope of its (complex) Fourier spectrum. Then we convolve its 
inverse Fourier transform with a white noise signal, in other words we create a stochastic signal with the same color as 
the irregular signal. Finally we mix this colored stochastic signal with the originally subtracted periodic component 
of the Rossler signal. The resultant stochastic data set is displayed in Fig. 1. Note that, by construction, the envelope 
of the Fourier spectrum of the stochastic signal is the same as that of the Rossler signal. 




100 200 ,. 300 400 

time 



FIG. 1. Input data sets: top: Rossler data; middle: stochastic data, Reconstructed signal: bottom: stochastic data 



By eye it is impossible to distinguish between the Rossler and the stochastic data sets. The Fourier spectra do not 
distinguish either because, by construction, they are very similar. What tools do we have to distinguish between the 
two signals? 

The study of many systems has shown that it is often revealing to make a 'first return map', i.e. to plot one 
extremum of the signal versus the following one . Fig. 2 compares the first return map for the Rossler (R) and 

the stochastic (S) signals. Clearly the S signal does not have the coherence of the R signal. The close resemblance of 
the R map to the logistic map ||^ suggests that the Rossler attractor is born out of a cascade of period doublings and 
that a horseshoe structure underlies its dynamics. We caution though that only very dissipative attractors can give 
rise to such nice compact first return maps. 
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FIG. 2. First return maps: left: Rossler data; right: stochastic data 
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While it is of course possible to plot the trajectory in a 3-D picture, this is generally not the optimal type of 
representation because often the attractor is compressed in some directions. It is preferable to project the trajectory 
onto Broomhead and King coordinates {^i}- These projections onto the eigenvectors of the correlation matrix 
have the property of providing an optimal spreading of the attractor in all 'directions'. The first column of Fig. 3 
shows the BK projections onto the first three axes for the R signal and the fourth column that of the S signal. Here 
we see again that the stochastic signal is fuzzier (less coherent). 
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FIG. 3. Broomhead-King projections; top: ^2 vs ^2, middle: ^1 vs ^2, bottom: ^2 vs ^3; Rossler data, col. 1; Reconstructed 
signal col. 2 3-D, col. 3 4-D; Stochastic data, col. 4, Reconstructed signal, col. 5 



Abarbanel et al. |g have developed a simple, yet very useful and powerful tool for estimating the embedding 
dimension, namely the method of 'false nearest neighbors' [Q. The underlying principle is best described by the 
following analogy: If we defined the distances between stars as being those which we measure on the celestial sphere 
we would find many nearest neighbors that are 'false' because they merely fall on neighboring lines of sight. We 
have to go to higher dimension (3 here) to correctly determine the distances. The method of false nearest neighbors 
generalizes this notion. One computes the distances between points X„ in higher and higher embedding dimensions 
ds- When the percentage of nearest neighbors does not change any more with increasing we have found the 
minimum embedding dimension. 

Fig. 4 displays the percentage of false nearest neighbors as a function of and the tolerance ( |^), and for two 
values of the delay parameter A. The percentages drop off very fast for the R set that is displayed in the two figures 
on the left, and they indicate that the R data are embeddable in 3-D. The situation is quite different for the S data 
set shown in cols. 3 and 4, indicating a much larger dimension (as it should because the dimension of a stochastic 
signal is indeed very large; the percentage still decreases with because for the small number of data points the 
density of points decreases so fast with dimension |^). 




FIG. 4. Percentage of false nearest neighbors: left: R signal, right, S signal 

A totally different characterization of an attractor, of a topological nature, has been discussed in Letellier et al. 
and, as an example, has been applied to the Rossler band as well. The idea is that a chaotic attractor can be 
described by the population of unstable periodic orbits, their related symbolic dynamics and their linking numbers. 
In three dimensional cases, periodic orbits may be viewed as knots and, consequently, they are robust with respect to 
smooth parameter changes allowing the definition of topological invariants under isotopy (continuous deformation). 
This topological approach is thus based on the organization of periodic orbits and symbolic dynamics. It is a very 



5 



powerful tool, but unfortunately, it only works for attractors for which the first return map consists of relatively well 
defined branches. 



B. The Flow Reconstruction 

We are ready now to discuss the global polynomial flow reconstruction. Our reconstruction contains several unknown 
parameters: First we have the dimension cIe- It is of course the minimal value of this parameter that we would like to 
determine, and we increase cIe until we obtain an embedding. Next we have the delay A. Studies of nonlinear time- 
series have shown that for many purposes there is an optimal time-delay A that is obtained from mutual information 
considerations (e.g. Q Q). However, for our analysis this value is not necessarily the best. For example, too short a 
delay collapses the attractor onto the diagonal and noise becomes a serious problem, whereas too large a delay causes 
a very strong nonlinearity in the map, such that a polynomial assumption cannot handle it. We therefore treat A 
as a free parameter in our search and anticipate (hope) that there is some range in which a polynomial map works 
fine. Finally, the parameter p is the highest order polynomial that we allow in the map. Again since the number of 
coefficients of the map goes like we would like to keep it to a minimum. 




FIG. 5. Error norm as a function of embedding dimension ds and polynomial order, p=2, ... 6, increasing from top to 
bottom 



In Fig. 5 we display a typical behavior of the error norm as a f unction of cIe and of p for fixed delay A=ll. We only 
show the error for the first component a—1 of X (c./. eq. ^.lOj ) because the other components behave similarly. As 
one would expect the error falls off monotonically with polynomial order p. In 3-D the error norm decreases relatively 
slowly though. A probable reason is that a low order polynomial approximation to the Rossler map may not be 
accurate enough (A rational polynomial map might be more appropriate, Q [pO{). 

The important result from Fig. 5 is that the norm shows a very clear drop and levelling off at dg around 3 and 4. 
This suggests that we should be able to construct a global map in 4-D, but that it may already be possible to do so 
in 3-D. This evidence alone is however not sufficient to establish the minimum embedding dimension, and we need 
to look also at the behavior of the properties of the synthetic traj ecto ries. By this we mean trajectories {X"} that, 
once we have constructed the map F, we obtain by iteration of eq. 2.4 starting with some seed value X". We refer to 
any component of these trajectories as a synthetic signal. Indeed it might be that the error norm is very small, 
indicating a good fit, but that the generic trajectories of the map are very different from the test data. 

Now a word about Lyaponov exponents (cf. e.g. |^) which give a useful quantitative measure of the map and of its 
attractor. In the case of the Rossler system where we know the 3 ODEs we can compute the 3 Lyapunov exponents 
in our chosen (a, b, c) regime directly - They are found to have the values Ai=0.057, A2=0, A3=-4.5, and 2.013. 
One notes that one of the Lyapunov exponents of the Rossler system is exactly zero because the system of ODEs is 
autonomous (a time independent flow) Next, the presence of a positive Lyapunov exponent indicates that the 
attractor is chaotic. Finally, the fractal Lyapunov dimension implies that it is a strange attractor 

The question now arises whether it is possible to recover the values of these Lyapunov exponents from the R test 
data? The answer is two-fold: First, the computation of all of the Lyapunov exponents from a time series necessitates 
long data sets, often at least an order of magnitude longer than is available from observations. Here, for example, the 
length of our R test-data set (4 000) is far too small to yield exponents that come even close to the correct values. 
In observational data, in addition, noise easily vitiates the reliability of the Lyapunov calculations. Second, however. 
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we can first construct the map and then generate synthetic trajectories of sufficient length to be able to compute the 
Lyapunov exponents. The reason why, with this two step procedure, we are not 'getting something for nothing' is 
discussed in ref. pH . 

Table 1 compares these Lyapunov exponents and Lyapunov dimension for synthetic trajectories for various values 
of embedding dimension and time-delay. In de>A the largest Lyapunov exponents (Ai) are seen to be very close to 
the exact value. An exception occurs for d^—i where Ai and A3 are too large by a factor of about 2. This fact is 
probably connected with the relative lack of stability of the map in 3-D. In the next section we shall see that the 
addition of small amounts of noise stabilizes the map and produces correct Lyapunov exponents. The scatter of the 
negative exponents is higher than for the positive ones, but this is to be expected 0. Table 1 shows that, except for 
dE= 3, the Lyapunov dimensions are close to their known exact values. For the sampling rate that we use here the 
exponents are essentially the same whether we compute them for the flow or for the map. 

We draw attention to the fact that d^ is independent of the embedding- dimension de- If our attractor reconstruction 
is to make sense this is of course a necessary condition, and conversely it is a check on our calculation. 

The second Lyapunov exponent A2 is always very small (« 10^"*). Although we have constructed a map here, 
the sampling-time is sufflciently small so that the map closely approximates a flow. The s mall ness of one Lyapunov 
exponent constitutes a very powerful check on our a priori assumption that a flow(eq. 2.4 indeed underlies the 
dynamics. 

As expected there is a range of delays for which one obtains good maps. In Fig. 1, bottom, we show the best 
synthetic signal (obtained with d^; = 4, p = 4). Is has very similar properties to the original set. (Because the 
time-series are chaotic, it is not very meaningful to compare them directly. Since the signals are so short, statistical 
comparison methods are not very useful here either. We therefore have to rely on a number of comparison tests, 
neither of which is totally convincing by itself, but when taken together they give a good level of confldence |l^] 
p]| . We note here that the Fourier spectrum of the synthetic signal is very similar to that of the R data set. The 
Broomhead King projections (Fig. 3) are also very similar. 

In contrast, the flow reconstruction of the S test data set yields only limit cycles, for all values of A and p. In fact 
the limit cycles are almost identical to the periodic part of the signal. The global flow reconstruction thus recognizes 
the stochastic noise for what it is and it does not erroneously indicate low dimensional behavior. 



C. Noise 



The reconstruction of a stable map encounters some difficulties for the Rossler data set. The reason behind it is 
that the neighborhood of the attractor is not well described. One can understand this as follows: Suppose we tried 
to apply this method to system which had a limit cycle. Unless we had a transient we would have no information 
whatsoever on the neighborhood of the limit cycle, and the neighborhood could equally well come out repulsive or 
attractive. Noise broadens the attractor and forces the map to form a stable neighborhood. 



TABLE 1. Lyapunov exponents and dimension, left: R set, right: R set + 0.2% gaussian noise 



ds A Ai A2 A3 A4 A5 dL \ ds A Ai A2 A3 A4 dL 



exact values: 



0.057 -4.5 2.013 



computed from synthetic signals: 



3 


5 


0.128 


<10-'* -6.4 






2.020 


3 


4 


0.046 


<10-'* -17.8 




2.003 


4 


4 


0.057 


-3.9 


-16.7 




2.014 


3 


5 


0.041 


-16.9 




2.003 


4 


5 


0.058 


-3.7 


-15.0 




2.014 


3 


6 


0.055 


-16.2 




2.003 


4 


6 


0.051 


-3.5 


-12.7 




2.014 


3 


7 


0.062 


-14.3 




2.004 


4 


7 


0.068 


-3.2 


-9.0 




2.021 


4 


4 


0.063 


-17.3 


-30.5 


2.004 


5 


4 


0.054 


-8.0 


-19.0 


-28.0 


2.007 


4 


5 


0.061 


-16.8 


-28.4 


2.004 


5 


5 


0.054 


-6.8 


-18.2 


-28.0 


2.008 


4 


6 


0.066 


-16.8 


-28.5 


2.004 


5 


6 


0.055 


-5.3 


-16.8 


-29.0 


2.010 


4 


7 


0.068 


-15.4 


-26.4 


2.005 


5 


7 


0.055 


-4.8 


-15.8 


-29.0 


2.011 
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We conclude from this study (1) that the minimum embedding dimension is 3 and (2) that K = int(di). From 
the inequahty ( ^.(j ) it thus foUows that d = 3. In other words the global flow reconstruction method recovers the 
dimension of the dynamical system ( |3.3| ) that has generated it. 



IV. A REAL SYSTEM: THE LIGHT CURVE OF R SCUTI 



In 1986 numerical hydrodynamic simulations of metal poor Cephcid models, called W Virginis stars, showed 
irregular pulsations ||2^ |Q . A further investigation revealed period doublings as a stellar parameter was varied 
and concluded that the irregular behavior of these models were a manifestation of low dimensional chaos. These 
calculations were based on models that ignored convective energy transport and treated radiation transfer in an 
equilibrium diffusion approximation. In order to confirm the theoretical prediction as to the chaotic nature of the 
irregular variability it was obviously necessary to compare to observations. 

Obscrvationally, it has been known for a long time that the light-curves of W Vir stars are regular (periodic) for 
low periods, but then develop irregularities as the period and luminosity increase. This increasingly irregular behavior 
carries over to their more luminous siblings, the RV Tauri stars. Unfortunately, these irregular variables have not 
received much attention from professional astronomers - the publication of what appears to be a table of random 
numbers is not compatible with career and funding pressures! Hopefully, the development of nonlinear time-series 
analyses such as this one will change this perception. However the Association of Variable Star Observers (AAVSO) 
in Cambridge collects observations and compiles data from amateur astronomers around the globe. Often these data 
are visual, and of poor quality. Yet, the large number of independent measurements make this a wonderful data base. 
For the variable star R Scuti a large data set of 15 years of data are available. The data have a normal distribution 
that is independent of magnitude, and it is therefore possible to extract a good light-curve despite the large individual 
errors. We show a short section of the observational data of R Set in Fig. 6 together with a smooth fit. It is this fit 
with equally spaced data points that we use as the input data for our global flow reconstruction. For more details of 
the smoothing procedure we refer the reader to the original papers [0 . The whole R Set light-curve is displayed 
in Fig. 7 on top. 




3000 3200 3400 3600 3800 4000 4200 4400 4600 4800 5000 

time [d] JD-2,440.000 
FIG. 6. Typical section of the observational data together with a smooth fit 



The Fourier spectral analysis of the light-curve shows a broad dominant peak at / ~ 0.007 d~^ and a secondary 
much broader structure centered on / w 0.014 d^"'^. The analysis of patchy, but very old data going back some 150 
years has shown that while the envelope of the Fourier spectrum is more or less steady, the individual peaks are not 
at all, especially the tertiary ones . There is a further discussion in ref . as to why the light-curve of R Scuti 
cannot be generated by multiperiodic pulsations, even when evolutionary changes are taken into account. 

Space does not permit us to describe the details of the global flow reconstruction and we refer the reader to the 
original papers [^| [H] (see also ||2^]). 

In Fig. 7, middle and bottom, we exhibit two pieces of synthetic signal (obtained mds =4 withp = 4). Visually, the 
synthetic time-series is very similar to the observational data. The Fourier spectrum similarly has the same envelope. 
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synthetic signal - cont 
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FIG. 7. R Scuti, top: Whole light-curve data set used in the analysis, middle and bottom: pieces of synthetic light-curves 
generated with the global map 



In Fig. 8 we show the Broomhead-King projections on the lowest three axes. The left-most colu mn represents the 
observational data. Columns 2 and 3 sho w ou r best reconstruction with 4-D and 5-D maps (eq. 2A). Finally, in col. 4 



we show a reconstruction with a flow (eq. 2.3) 



In Table 2 we show the Lyapunov exponents and dimension for several values of A and cIe- Noteworthy are the 
following points: (1) the results show a certain robustness with respect to the delay A, (2) ~ 3.1 is essentially 
independent of d^, (3) Ai > which establishes the chaotic nature of the attractor, and (4) A2 ~ 0, which confirms 
the presence of a flow. We note that it has been found impossible to construct a map with synthetic signals that 
bear any resemblance to the data which indicates that the minimum embedding dimension must be greater than 3. 
However, as Table 2 shows in 4-D we can construct a flow whose propeties remain largely invariant as we increase 
ds- We conclude that the minimum is therefore 4. We add in passing that a false nearest neighbor analysis 
corroborates 4 as the minimum embedding dimension. 



R Set obs Syn 4D Syn 5D Syn 4D F 




FIG. 8. Broomhead-King projections; top: ^2 vs ^2, middle: vs ^2, bottom: ^2 vs ^3; left: R Scuti observational data, 
cols. 2 and 3: map reconstructions in 4-D and 5-D, right: flow reconstruction in 4-D; 



Having now established from the Lyapunov dimension that the lower bound K+1 (rel. |^) of the physical phase- 
space dimension is 4 we thus conclude that d = 4, i.e. the light-curve of R Scuti is generated by a 4-D dynamics. 
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TABLE 2. Lyapunov exponents and dimension, R Scuti 





A 


P 


Ai 


A2 


As 


A4 


dL 


4 


4 


4 


0.0019 


<io-4 


-0.0016 


-0.0061 


3.05 


4 


5 


4 


0.0017 




-0.0014 


-0.0054 


3.06 


4 


6 


4 


0.0019 




-0.0009 


-0.0051 


3.19 


4 


7 


4 


0.0020 




-0.0011 


-0.0052 


3.18 


4 


8 


4 


0.0014 




-0.0010 


-0.0049 


3.07 


5 


7 


3 


0.0016 




-0.0005 
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V. CONCLUSION 



The polynomial global reconstruction method has been tested on the Rossler attractor. Using as input a short 
time-series of only one of the 3 Rossler variables one obtains a very good estimation of the Lyapunov exponents and 
the fractal dimension of the attractor. Furthermore it is possible to determine that this scalar signal is generated by 
a 3-D flow (the Rossler equations 3.3), in other words one can find out the dimension of the physical phase-space of 
the dynamics. 

Encouraged by these results the flow reconstruction has been applied to real data, viz the light-curve of the variable 
star R Scuti. Not only are these data contaminated with a good amount of observational and other extrinsic noise, 
but we have no a priori information about the dynamics that has generated them. 

The analysis has shown that the dynamics which generates the light-curve of R Scuti is of dimension four, i.e. the 



signal can be generated by a flow in 4-D (eq. 2.1). Furthermore the attractor has a fractal dimension « 3.1. This 



must be considered remarkable considering that this star undergoes luminosity fluctuations of up to factors of 40 (cf 
Fig. 1), with shock waves and ionization fronts criss-crossing the stellar envelope. Some of the physical implications 
have been addressed in references ||ll[] |^ . 

If nothing else it is interesting that the nature of the irregular variability of these types of stars is no longer a 
mystery, but that it has such a simple explanation in terms (very) low dimensional chaos. 
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